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Abstract 

The electrical conductivity of graphene containing point defects is studied within the binary alloy 
model in its dependence on the Fermi level position at the zero temperature. It is found that the 
minimal conductivity value does not have a universal character and corresponds to the impurity 
resonance energy rather than to the Dirac point position in the spectrum. The substantial asym- 
metry of the resulting dependence of the conductivity on the gate voltage magnitude is attributed 
as well to the very shift of the conductivity minimum to the resonance state energy. 

PACS numbers: 71.23.-k, 71.55.-i, 81.05.ue 
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I. INTRODUCTION 



Graphene, a thermodynamically stable graphite monolayer, which has been mechanically 
exfoliated for the first time only a few years ago,— - — is gaining considerable scientific attention. 
This new material looks promising enough for a number of important practical applications, 
some of which have been long hoped for. While experimenters are targeted at engineering 
graphene based devices in the not so distant future, graphene attracts theoreticians as the 
first existing in the free state physical system, which can be named two-dimensional (2D) 
without any reservations. Undoubtedly, so far unique electronic properties of graphene were 
the most challenging issue. These properties directly come up from the honeycomb lattice 
with its two-atomic structure, inherent in a single atomic layer of graphite. The lattice 
structure leads to the Dirac dispersion of charge carriers, which makes up the core of studies 
devoted to graphene. 

Transport properties of this material are, sure enough, of primary importance for 
graphene-based electronics. In real crystals, transport properties essentially depend on 
non-ideality of the system and on interaction of carriers with other excitations. Below we 
are going to focus on imperfections of graphene, and, in particular, on point defects in it. 
This allows to rise a question on the spectrum of delocalized carriers, on its dependence on 
the amount of defects, and, eventually, on such a remarkable quantity as the minimal value 
of the conductivity. 

Initially, two main features of the graphene conductivity were singled out. First of all, the 
conductivity of graphene devices never dropped below a certain value. Since this minimal 
value seemed not to vary between different experimental samples, the origin of the universal 
behavior of the minimal conductivity value has been extensively searched for. These efforts 
shaped the famous minimum conductivity puzzle. The linear dependence of the conductivity 
on the gate voltage made up the second feature. However, the minimal conductivity value has 
been found soon to be strongly sample dependent,- and the effect of minimum conductivity 
has been attributed to graphene's imperfections. In view of mentioned features of the 
graphene conductivity, a qualitative difference between charged impurities and point defects 
has been established.- While it has been demonstrated that charged impurities are able 
to yield the required linear dependence of conductivity on the gate voltage, point defects 
were shown to produce the sub-linear conductivity behavior, and, consequently, ruled out as 
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the conductivity limiting factor. The concept of charged impurities as a main source of the 
scattering of charge carriers in graphene has been thoroughly developed,- - - and convincingly 
compared to the experimental data on graphene with deposited potassium atoms on its 
surface.-* 1 ^ At that, just a constant contribution to the conductivity were ascribed to point 
defects when fitting the experimental data. 

Even though the concept of charged impurities looks sounding, experiments on graphene 
in ethanol environment seriously question the dominant role of the Coulomb scatterers.— In 
addition, the conductivity asymmetry evident in measurements of graphene with deposited 
potassium atoms has not received the proper explanation yet.-* 1 ^ While a moderate asym- 
metry can be attributed to the disbalance of positively and negatively charged impurities,— 
the marked asymmetry of conductivity dependence on the gate voltage in graphene doped 
by transition metals manifests the response that is different from the one, which is expected 
from the charged impurity centers.— The clearly sub-linear character of conductivity curves 
corresponding to graphene samples heavily doped by transition metals or lightly doped by 
the potassium atoms only strengthen the overall impression that we are dealing with the 
interplay of different types of disorder, and that each one of them should receive a com- 
prehensive treatment in an effort to grasp the conductivity properties in graphene. As an 
example, an interplay between charged impurities and point defects, which involves differ- 
ent reaction of these two impurity types to screening effects, can be employed to solve the 
dilemma of slightly varying minimal conductivity value for graphene placed into a dielectric 
environment. 14 

It must be stressed that in the mentioned experiments targeted at finding the main 
scattering channel for carriers in graphene, different adatoms were deposited on the surface 
of samples, which otherwise were considered as pristine. In fact, here we are dealing with 
two different issues again. One of them is determining what actually limits the conductivity 
of graphene samples obtained by a certain technique, and another one is analyzing the effect 
of intentionally added impurities on the conductivity features. The latter is closely related 
to the current tendency to functionalize graphene by substitutionals, adatoms, or chemically 
active groups. Regarding the conducting properties, such an adjustment of graphene can 
proceed up to a possibility of the metal-insulator transition, which has been successfully 
observed in graphene doped by hydrogen atoms recently— Similarly, graphene demonstrated 
insulating behavior after irradiation by Ne ions, which is expected to produce short-range 
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defects.— 

In view of that, before any complex models of impurity centers are constructed, the simple 
ones should be properly examined in their characteristic aspects. And the basic model of 
the point defect is definitely among them, since it allows for a Mott transition in impure 
graphene.— 

Below we are returning to the common in semiconductor physics model of a binary alloy 
intending to examine what features of the graphene conductivity it is capable to reproduce. 
Such short-range impurities violate the electron-hole symmetry of the system and are nat- 
urally providing for the conductivity asymmetry. This impurity model has been studied 
either in a weak scattering limit, or in the unitary limit. We are going to show that this 
model does exhibit some worthwhile features in-between these two extremes. In particu- 
lar, we demonstrate that the effective conductivity minimum, in contrast to former studies, 
corresponds not to the Dirac point of the spectrum, but to the energy of a single impurity 
resonance, where the impurity scattering is the strongest. 



II. MODEL DISORDERED SYSTEM 



The host Hamiltonian for electrons in graphene is taken in the tight binding approxima- 
tion with hopping restricted to the nearest neighbors, 

H = t Yl c n«c m/3 , (1) 

<na,m/3> 

where t f» 2.7eV is the hopping parameter for the nearest neighbors,-^ n and m run over 
lattice cells, a and (5 enumerate two sublattices of the honeycomb atomic arrangement, $ na 
and c na are the electron creation and annihilation operators at the respective lattice cites. 
Substitutional impurities are supposed to be distributed evenly and in uncorrelated manner 
on lattice sites. Presence of an impurity at a given lattice site is assumed to be manifested 
only through a change in the respective on-site potential of the tight-binding Hamiltonian. 
This type of impurity perturbation fully corresponds to the conventional model of a binary 
alloy with a diagonal disorder, which had been extensively used in physics of real crystals, 
and sometimes is referred to as the Lifshitz (or isotopic) model for historical reasons. The 
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corresponding Hamiltonian for a disordered graphene has the form, 

H = H + H imp , H imp = V L ^2'i] na cl ia c na , (2) 

where Vl is the deviation of the potential at the impurity site, and variable r\ na is unity 
with the probability c or zero with the probability (1 — c), which specifies c as the impu- 
rity concentration. The "per carbon atom" concentration c can be easily converted to the 
impurity coverage, 

n im = n c, (3) 

where 

is the inverse area per one carbon atom, and a = 0.246 nm is the lattice constant of graphene. 
The Green's function of a disordered system, 

g = (E- H)-\ (5) 

after averaging over different impurity distributions, G = (£/), regains the translational 
invariance and can be expressed by means of the Dyson equation, 

G = <7 + </£G, (6) 

where £ is the self-energy, and 

g(E) = {E- Ho)' 1 (7) 

is the host Green's function. When the amount of introduced impurities is moderate, it is 
possible to implement the modified propagator method.— Within this approach, the self- 
energy is site diagonal and identical on both sublattices, 

E.W, ^ BW ^ = __£^_ M , (8) 

where J is the identity matrix. At a small impurity concentration, c C 1, multiple occupancy 
corrections are not significant. Thus, the method of modified propagator yields results 
that are practically indistinguishable from the ones produced by the conventional coherent 
potential approximation. 
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The calculation of the diagonal element of the host Green's function in the site represen- 
tation g na na, which will be required for the subsequent analytical treatment of the impurity 
problem, is given in Appendix. From here and on we are choosing the bandwidth parameter 
W, see Eq. 0A.9p . as the energy unit. For the sake of clarity, we will use for the dimensionless 
energy and the impurity potential following designations: 

< = #• » = £ < 9 > 

Since only energies small compared to the bandwidth, i.e. those, at which the linear dis- 
persion holds in the host system, are considered, the corresponding approximation to the 
dimensionless diagonal element of the Green's function can be written as follows, 

Wg nana (E) =g {e) « 2eln|e| -i7r|e| , |e| < 1. (10) 

It is not difficult to see that the diagonal element of the Green's function f JTOj) looks 
similar to the properly scaled diagonal element obtained within the model of massless Dirac 
fermions for the electron spectrum in graphene. Thus, despite the fact, that the intervalley 
scattering is, for sure, taken into account in Eq. (JHJ) explicitly, the net result for the single- 
site scattering will not be qualitatively different from the one for the frequently used model 
of massless Dirac fermions, in which only one Dirac cone is retained. This resemblance 
follows from the single-site character of the impurity perturbation (see Eq. (j2J)). At that, 
the validity of the modified propagator method (jSJ) is limited by the scatterings on impurity 
clusters,— 1 ^ which contribution to the self-energy will be monitored in what follows. 



III. RENORMALIZED ENERGY PHASE AND CONDUCTIVITY 

In order to make the self-consistency condition (|SJ) more tractable, a regular substitution 
can be made, 

e — S(e) = xexp(^), x > 0, < if < ir, (11) 

which singles out the phase of the renormalized energy e — S(e). This phase diminishes 
from 7r/2 to zero inside the conduction band and rise from tt/2 to tc within the valence band 
when moving away from the Dirac point position. With the help of the obtained above 
expression ( TTOl for the diagonal element of the Green's function and the substitution ( flTi) . 
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the imaginary part of Eq. (jSJ) can be reduced as follows, 



cv 2 [2 In x + (2ip — 7r) cot ip] + [1 — t>x(2 In x cos ip — (2<^ — 7r) sin<^)] 2 + 

+ [vx(21nxsiny2 + (2<p - vr) cos^)] 2 = 0. (12) 

Provided that the impurity perturbation strength v and the impurity concentration c are 
fixed, this equation establishes a correspondence between the renormalized energy modulus 
x and its phase (p. For those x that are exceeding a certain threshold magnitude, which 
is, indeed, determined by the impurity concentration and the perturbation strength, this 
equation always has two different solutions with respect to the phase (p. One of them 
(ip < 7r/2) belongs to the conduction band, while the other (ip > n/2) lies within the 
valence band. The literal carrier energy that corresponds to a given renormalized energy is 
determined by the real part of Eq. (jSJ) , 



e = x cos ip+ 



+ 



cv [1 — vx(2 In xcos ip — (2<p — 7r) sin ip)] 



[1 — t>x(2hixcos<y2 — (2ip — 7r) siny?)] 2 + [vx(2\ia xsimp + (2<p — it) cosy?)] 2 



(13) 



Taken together, the last two equations, f lT2|) and (TIB"]) , are making up a set, which implicitly 
specifies the dependence of the renormalized energy phase ip on the carrier energy e. 

Since the procedure required to calculate the self-energy is already outlined, it is possible 
to employ the Kubo expression for the conductivity of a disordered graphene at the zero 



temperature 
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e F — ReS(e^)^ — ImS(e i ?) 



arctan 



e F — Re £(ep) 
— Im S(ep) 



(14) 



— ImS(e j p) ' 6p — ReS(e^) 

where e F is the Fermi energy. By means of the substitution (llip . which was used above to 
simplify the self-consistency condition for the self-energy (jSJ), the above expression can be 
significantly reduced, 

,2\ 2 



&cond \ 7 | ^condi 



"cond 



1 + (cot <^>i? + tan </?p) ^— — ipf^j 



(15) 



where is the renormalized energy phase at the Fermi level, and the dimensionless con- 
ductivity cr con rf, which will be used onwards, is singled out. It should be emphasized that 
the dimensionless conductivity u CO nd depends on the renormalized energy phase ip alone. In 
the same way, the well-known Ioffe-Regel criterion,— which is commonly used to separate 



extended states in a disordered system, and the applicability criterion of the modified prop- 
agator method can both be expressed through the same renormalized energy phase.— It has 
been shown that with varying the renormalized energy phase the modified propagator ap- 
proximation validity violation and the indication of the state localization by the Ioffe-Regel 
criterion are occurring simultaneously for those states, which energies fall inside the host 
band of a disordered system.—^ Certainly, it is not conceptually correct to expect that the 
Ioffe-Regel criterion can be used to pinpoint precisely the mobility edge position in a disor- 
dered system. Similarly, there should be no sharp boundaries between those spectral regions, 
in which the modified propagator method is applicable, and those ones, in which it is not. 
Nevertheless, there are strong arguments supporting the estimation that the mobility edge 
in a disordered system should be located at those energy, at which the renormalized energy 
phase is close to 7r/6 for the conduction band, and, respectively, to 57r/6 for the valence 
band. Thus, in those spectral intervals, inside which states are anticipated to be localized 
according to the Ioffe-Regel criterion, neither the Kubo formula ( Fl5l) has any relevance, nor 
the modified propagator method is reliable. On the contrary, the approach outlined above 
is consistent in the spectral domains occupied with extended states, where the renormalized 
energy phase p is either small (for the conduction band) or close to tt (for the valence one). 



IV. CONDUCTIVITY IN DIFFERENT SCATTERING REGIMES 



A. Weak scatterers 



When the impurity perturbation strength is moderate (\v\ < 1), it is possible to take an 
advantage of the renormalized energy phase smallness (or its closeness to tc) and construct 
a correspondent approximate solution of Eq. f fT2|) . 
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d ~ (1^2t;xlnx) 2 + (7n;x) 2 + 2ct; 2 (l + lnx)' ° < ^ 
where 9 stands for (p inside the conduction band, and for ir — tp inside the valence band. 
The sign in the denominator also switches from a minus to a plus when moving from the 
conduction band to the valence band. Obviously, the renormalized energy phase is close 
to 7r/2 in a narrow interval of energies around the shifted Dirac point, and thus the above 
approximation is not valid inside this region. However, the transition of the renormalized 
energy phase from small values to values that are close to it is very fast. This transitional 
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region is, in fact, exponentially narrow and, for certain reasons, should be treated separately, 
as it will be explained in detail below. 

It is not difficult to check that in this scattering regime (\v\ < 1) the effective shift of 
states along the energy axis, which is given by the real part of the self-energy ReS(e), is 
nearly constant in the whole domain under consideration (|e| <C 1). Therefore, as a first 
approximation, one can take 

±xsie-c?;, (17) 

where the sign is varying according to the current band as above, so that x always remains 
positive, as it should do. The expression for the conductivity, Eq. f lT5|) . can be also simplified 
utilizing the smallness (or closeness to 7r) of the renormalized energy phase, 

(Tcond « ~, < 1, (18) 

where linear terms and terms of the higher order in 9 are omitted. All these approximations, 
Eqs. (TIS]) - (118I) . can be combined into the final expression for the dimensionless conductivity, 



[1 - 2v (e F - cv) In \e F - cv \} 2 + \irv(e F - cv)} 2 2 ri n . . 

Vcond-- — — ' 2 — + - 1 + ln e F -cv , 19 

7CCV 1 71 

which fits well the conductivity, calculated numerically by Eqs. (fT!?j) .f fT3"j) . and ffTBT) . through- 
out the whole considered interval of energies (|e| <C 1). 

As follows from (fT9j) . the conductivity is gradually diminishing with increasing the Fermi 
energy e F from the valence band to the conduction band for a negative impurity perturbation 
(v < 0) and vice versa. The conductivity of graphene calculated by Eqs. (jT2l) . ([13]) . and (|T5|) 
without any additional approximations at different concentrations of point defects is plotted 
against the Fermi energy in Fig. [T]for the case of moderate impurity perturbation (\v\ < 1). 
On the whole, the dependence of the conductivity on the Fermi energy is smooth and almost 
featureless, while being strongly asymmetric against the shifted Dirac point position. The 
only exception from the monotonic behavior of the conductivity can be observed in a close 
vicinity of the Dirac point, where the curve manifests a sharp dip, which is barely discernible 
in the curves corresponding to high impurity concentrations. If to trust the results yielded 
by the modified propagator method all the way down to the Dirac point, the conductivity 
at its very tip should drop to 4/rr. This directly corresponds to the notorious theoretical 
magnitude of the universal minimum conductivity in the inhomogeneous graphene, which 
has been widely debated within the so-called "missed 7r" discourse. However, the modified 
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FIG. 1: Conductivity of graphene with point defects vs Fermi energy for v = —0.5 and concentra- 
tions c = 0.1/2™, 7i = 1 ... 5. Arrows point at positions of narrow dips. 

propagator method is not applicable in the Dirac point neighborhood. It can be shown^ 1 ^ 
that this approximation is not reliable in the interval 



Therefore, as it was outlined in the previous section, the Kubo formula is also ineffective in 
this interval, and the obtained conductivity magnitude at the Dirac point has no physical 
meaning. Still, the width of the energy interval, in which the analytical approach fails, 
is exponentially small compared to the bandwidth. Since the corresponding dip on the 
conductivity curve is so narrow, it should be averaged out at realistic sample temperature, 
or by means of any other broadening mechanism. 

Consequently, the presence of the sharp dip on the conductivity curve can be neglected, 
and, probably, should never come out in actual experiments. The asymmetry of the con- 
ductivity dependence on the Fermi energy arises from the presence of the severely smeared 
out impurity resonance, which enhances the impurity scattering. Indeed, the smooth char- 
acter of the conductivity curve does not resemble the experimentally observed check mark 




(20) 
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shape. However, if this check mark shape is caused by another dominating type of impu- 
rities, weakly scattering point defects undoubtedly can contribute to the asymmetry of the 
conductivity curve. 



B. Strong scatterers 

In the limit of the strong impurity scattering, \v\ 3> 1, situation is completely different. 
When the impurity potential is large compared to the bandwidth, a well-defined resonance 
state is manifested in the electron spectrum.— In the limit of a strong impurity potential, 
the resonance state energy e r is determined by the Lifshitz equation, 

1 ?y v Reg (e r ) ps 2ve r In \e r \ , (21) 

while the resonance state damping is given by 

r r « 7r| f rl (22) 



For the resonance state to be well-defined, the condition 

* = rr ~ on . i i I, « 1 23 

\e r \ 2 |1 + In |e r || 

must be met. Thus, one should have | In |e r || ^> 1, which corresponds to a strong impurity 
perturbation and a resonance energy located close to the Dirac point. 

The qualitative difference of the strong impurity perturbation case resides not only on 
the presence of a resonance state in the spectrum, but mainly on the fact that in this case 
the electron spectrum undergoes a radical rearrangement. That is, with increasing impurity 
concentration a quasigap filled with localized states opens up around the resonance state 
energy . — i i & M There exists a certain critical concentration of impurities, 

which is determined by the mutual spatial overlap of individual impurity states. When 
impurity concentration exceeds the critical concentration c r of the spectrum rearrangement, 
the width of the quasigap starts to increase rapidly with increasing impurity concentration 
as y/ — 2c/ In c— Certainly, neither the modified propagator method nor the Kubo ex- 
pression for the conductivity will work inside this quasigap. Therefore, we will consider only 
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FIG. 2: Conductivity of graphene with point defects vs Fermi energy for v = — 2 and concentrations 
c = c r /2 n , n = 1 . . . 5, c r ~ 0.012. Arrows point at positions of narrow dips. 

those impurity concentrations that are less than the critical one (c < c r ) in the case of the 
strong impurity potential. Vacancies are frequently modeled by point defects with infinite 
impurity potentials v. Because of this, the critical concentration c r for vacancies in graphene 
is zero. In other words, the spectrum rearrangement is already over for any concentration 
of vacancies. Therefore, vacancies are out of the scope of the present study. 

The conductivity calculated directly by Eqs. ( JT2l) . ( TL3l) . and ( fi~5l) at different concen- 
trations of point defects is shown in Figs. [2] and |3] for a not so excessive (v = —2) and 
for a reasonably strong (v = —8) impurity potential, respectively. The Dirac point shift, 
which occurs along with the impurity concentration increase, is not so pronounced, because 
c r v ~ 1/v. Like in the case of the weak impurity potential, there is a sharp dip in the 
conductivity curve located at the Dirac point. The hint of this dip can be seen in the figures 
at concentrations that are approaching the critical one. Nevertheless, the presence of this 
dip should be neglected by the same arguments as above. 

What really distinguishes the strong impurity perturbation case is the presence of the clear 
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FIG. 3: Conductivity of graphene with point defects vs Fermi energy for v = — 8 and concentrations 
c = c r /2 n , n = 1 . . . 5, Cr « 0.0005. Arrows point at positions of narrow dips. 

minimum on the conductivity curve, which is located at the energy of the impurity resonance 
state. This is understandable, since the impurity scattering is the strongest around the 
resonance energy. The width of this minimum corresponds to the resonance state broadening, 
and, therefore, this minimum is not as sharp as the minimum at the Dirac point. Overall, the 
conductivity curve acquires a quasi-parabolic form, which is particularly well-pronounced at 
lower impurity concentrations. In addition, the conductivity curve appears more symmetric 
for a larger impurity potential. 

With increasing the impurity concentration, the concentration broadening of the reso- 
nance state also increases. At the impurity concentrations that are close to the critical 
one, the broadening of the resonance state is as wide as the distance from the resonance 
energy to the Dirac point. This widening of the resonance broadening area along with the 
tendency of states toward localization inside it are manifested by the apparent flattening of 
the conductivity curve around the resonance energy at c ~ c r . Outside the domain of the 
concentration broadening, the approximate expression for the conductivity is even simpler 
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than before, 



[1 -2w F ln|e F |] 2 + [irve F ] 2 

O'cond ~ 5 • (25) 



It is not difficult to see from Eq. ( !25l) that in the unitary limit 



acU _ „ (26) 

which corresponds to the known result for vacancies in graphene.^ As was stated above a 
quasigap around the Dirac point should be present at any impurity concentration in the 
unitary limit of the impurity perturbation. The approximation f[2T)j) is valid only outside of 
this quasigap. 

For the finite impurity perturbation, the expression ( 125]) also can not be used close to 
the resonance energy. In order to obtain the minimum value of the conductivity, it is 
required to know the magnitude of the renormalized energy phase at the resonance energy. 
Its concentration dependence, (p r {c), follows from the self-consistency condition ( fl2|) . The 
second term in this equation nullifies by the very definition of the resonance energy Eq. (1211) . 
The remaining two terms constitute the relation: 

c— — 2e^(c) tan<yj r (c) In — - — — taxLcpJc) + <£v(c) — — ■ (27) 

L cos (p r (c) 2- 

This expression can be significantly simplified by taking into account that introduced earlier 
symmetric phase 9 is small at low impurity concentrations, 

c^7ie 2 r (c)6 r {c)\l + ^-], P (c)<l. (28) 
L 7 r J 

Dependence of the resonance energy on the concentration e r (c) is very weak for the strong 
impurity perturbation and can be neglected. By setting the resonance energy in Eq. ( 128]) to 
its value for the isolated impurity, this equation can be solved for the renormalized energy 
phase at the resonance 9 r , 



w«! V^? -1 )• r —^r (29) 



,2^ 



The minimum conductivity at the resonance energy is then given by Eq. ( 1181) . Indeed, 
it does not again have the universal character and varies with impurity concentration. The 
minimum value of conductivity calculated numerically with the help of Eqs. (fT2l) . (fT3|) and 
([15]) is plotted against the impurity concentration in Fig. H] for two different values of the 
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FIG. 4: The minimum conductivity vs impurity concentration for v = —2 (upper curve) and for 
v = —8 (lower curve). Corresponding analytical approximations are given in dashed lines. 

impurity potential. Initial fast minimum conductivity drop, which occurs with increas- 
ing the impurity concentration, is followed by a considerable flattening of the curve. The 
manifested saturation-type behavior of the minimum conductivity concentration dynamics 
qualitatively corresponds to the observed data.- According to Eqs. ( 129|) and ( !T8|) . at small 
impurity concentrations, c<c*, the minimum conductivity of graphene with point defects 
is proportional to 1/c, which is similar to the case of charged impurities.^ However, at 
higher impurity concentrations, c 3> c*, the minimum conductivity for point defects falls 
down more slowly, namely as 1/ i/c. Thus, if both expressions for the minimum conductivity 
are fitted to each other in the low concentration limit, then the one corresponding to the 
point defects should yield significantly larger values of the minimum conductivity at c ^> c*. 

V. CONDUCTIVITY ASYMMETRY 

The conductivity of graphene devices is usually measured against the applied gate voltage. 
Since the gate voltage controls the carrier density in the graphene sample, these experimental 
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curves can be simulated by plotting the conductivity as a function of the number of occupied 
states. Leaving out the irrelevant constant, and taking into account two actual sublattices, 
the number of occupied states can be written as follows, 

n (e F ) = -- [ F Im({e - S(e)}{2 ln[e - E(e)] - m})de. (30) 

Jo 

The conductivity of impure graphene, calculated as before by Eqs. ( |T2|) . (113"]) and ( !T5|) . is 
plotted in Figs. [5] and [6] for both chosen strengths of impurity perturbation against the 
number of occupied states, which is given by Eq. ( 13"U|) . The number of occupied states is 
calculated by the numerical integration. The change in the introduced magnitude An(e^) 
can be easily related to the respective change in the carrier density An, 

An = n An(e^), (31) 

where n is defined by Eq. (TJJ. In the usual experimental setup, the carrier density depends 
linearly on the gate voltage V 9 , so An = xv^V 9 , where xv ~ 7.3 x 10 10 cm~ 2 V -1 . Thus, 
when the gate voltage is varying in the window of ±100 V, the value of n{ep) varies by 
0.004. 

It is visible from Figs. [5] and |6] that the calculated conductivity dependence on the gate 
voltage is highly asymmetric. Similar asymmetric character of the conductivity curve has 
been already reported elsewhere for the graphene with point defects. 27 While the expected 
slightly sublinear behavior of the impure graphene conductivity is readily reproduced, the 
asymmetry of the curve appears to be a bit on the extreme side. Although such a strong 
asymmetry is sometimes reported for a graphene with deposited adatoms, 13 its origin for 
the point defects is to be understood. In order to proceed in this direction, the conductivity 
can be expanded into a series in the vicinity of its minimum. Since it has been reasoned 
above that the conductivity reaches its minimum value at the resonance energy e r for the 
strong impurity perturbation, the expansion is straightforward, 

a cond ^a° cond + Xr(e F -e r ) 2 } Xr > (32) 

where cx° ond is the minimum value of the conductivity, and Xr is some constant. 

In the current paper we restrict ourselves for the strong scattering regime to those im- 
purity concentrations that are less than the critical concentration c r of the spectrum rear- 
rangement. In this case, the density of states is not considerably distorted by the presence 
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FIG. 5: Conductivity of graphene with point defects vs number of carriers for v = —2 and concen- 
trations c = co/2 ra , n = 1 . . . 5, co ~ 0.012. 

of defects. Because of the estimative character of this arithmetic, it is quite sufficient to 
assume that the density of states remains completely unchanged, i.e. identical to the host 
system, except of the rigid shift of both bands to a new Dirac point e^, 



V 9 « V 9 D + X D(tF - t D f sgn(e F - e D ), X d > 0, 



(33) 



where Vp is those magnitude of the gate voltage V 9 , at which the Fermi level comes to the 
Dirac point of the spectrum, and xd is some constant. This equation can be easily solved 
for the Fermi energy, 

l\V9 - V£\ 

(34) 



|r ' V -^sgn(V 9 -V 9 ). 



Xd 



e F ^e D + 

Substituting this result to the expansion (1321) . one can obtain 

X: 



<?cond ~ ^cond + 



XD 



V 9 -V^sgniV 9 -V 9 ,) 



XM e r - £d) 



(35) 



It is evident from this relation, that the conductivity asymmetry arise from the shift ex- 
pressed by the second term in the square brackets. If the conductivity reaches its minimum 
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" 0.002 0.004 0.006 0.008 

n(e F ) 

FIG. 6: Conductivity of graphene with point defects vs number of carriers for v = —8 and concen- 
trations c = co/2 n , n = 1 . . . 5, co ~ 0.0005. 

precisely at the Dirac point, then the conductivity dependence on the Fermi energy is linear 
and symmetric. However, we have demonstrated that the conductivity minimum is attained 
at the impurity resonance energy, which is essentially different from the Dirac point energy. 
This very difference does form the ground for the substantial conductivity asymmetry. 

VI. CONCLUSION 

It is demonstrated that there are two scattering regimes, which characterize the behav- 
ior of the conductivity in graphene with point defects. In the weak scattering regime, i.e. 
when the impurity perturbation strength is less than the bandwidth, the dependence of the 
conductivity on the Fermi energy is monotonic and asymmetric, which can contribute to 
the observed conductivity asymmetry, when point defects does not dominate other sources 
of scattering. However, in the strong scattering regime, i.e. when the impurity potential 
exceeds the bandwidth, the conductivity caused by point defects manifests a distinctive 
minimum in its dependence on the Fermi energy. This minimum, in contrast to a majority 
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of anticipations, corresponds not to the Dirac point of the spectrum, but to the impurity 
resonance energy. In this regime, the asymmetry of the conductivity dependence on the 
Fermi energy is noticeable. What is more, the pronounced asymmetry of the corresponding 
dependence of the conductivity on the gate voltage is caused by the very shift of the conduc- 
tivity minimum from the Dirac point to the impurity resonance energy. Despite the basic 
nature of the considered impurity model, it can qualitatively capture the essential features 
of the impure graphene conductivity manifested in experiments. Thus, one can expect that 
increasing the number of parameters characterizing the point defect will permit to approach 
closer to the quantitative description of conductivity features in graphene with point defects. 
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Appendix: Diagonal element of Green's function 

The straightforward expression for the diagonal element of the host Green's function 
reads: 

i r e 

9nana (E) = ^-J E2 _ E2{k) dk, (A.1) 
where the integration is carried over the entire Brillouin zone, which has the area 

Sbz = ^ (A.2) 

and E(k) is the unperturbed dispersion relation corresponding to the host Hamiltonian (OQ). 

In practical situations the Fermi level in graphene is located, nearly unavoidably, in a 
narrow spectral region, in which the dispersion is linear with a good accuracy. Near each of 
the two inequivalent Dirac points, the dispersion relation E(k) can be expanded, 

E(k') « ±v F k', (A.3) 

where k' is taken relative to the corresponding Dirac point, and 

vf = (A.4) 
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is the Fermi velocity. Consequently, the integration in Eq. (lA.lj) over the wave vector can 
be also performed relative to the each Dirac cone vertex, 



9o(E) 



E 



-dk. 



E < 3t, 



(A.5) 



S BZ J E 2 -E 2 (k) 

where the factor of 2 reflects the existence of two Dirac cones in the spectrum. However, due 
to the mutual overlap between respective cones, the integration in ( 1A.5I) can not be done 
over the entire Brillouin zone within the linear approximation for the dispersion relation. 
The corresponding cutoff magnitude of the wave vector is determined by the sum rule, 



47T 



s 



BZ Jo 



kmax 



dk 



which yields 



2^F 



\/3c 



(A.6) 



(A.7) 



Then, the integration can be performed exactly, 



90(E) 



47T 



E 



Sbz Jo & ~ (v F k) 2 



kdk 



E 



o E 2 -V37it 2 x 

J2 



dx 



€ 

W 



In 



1-e 2 



m sgn e 



where 

is the bandwidth parameter for the pure Dirac spectrum. 



E <W, (A.8) 



(A.9) 
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